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The dynamical scaling of ageing ferromagnetic systems can be generalized to a local scale 
invariance. This yields a prediction for the causal two-time response function, which has 
been numerically confirmed in the Glauber-Ising model quenched into the ordered phase. For 
a dynamical exponent z = 2, a new embedding of the Schrodinger group into the conformal 
group and the resulting conditions for the validity of local scale invariance are discussed. 



1 Phenomenology of ageing in simple ferromagnets 

Ageing phenomena provide a paradigmatic example of collective behaviour far from equilibrium 
and have received a lot of attention in recent years IHl • Ageing has been observed first 
in glassy systems, but for an improved conceptual understanding, it might be more useful to 
study first ageing phenomena in the simpler ferromagnetic systems, as we shall do here. We 
consider a ferromagnet with a critical temperature Tc > and prepare the system in some initial 
state (which typically may be disordered). Then the system is suddenly brought into contact 
(quenched) with a heat bath of temperature T < Tc (or T < Tc). Keeping T fixed, the system's 
evolution towards its equilibrium state at temperature T is observed. For definiteness, consider 
an Ising-like system with a microscopic degree of freedom cr^ = ±1 which can only take two 
possible values, and where i denotes the sites of a (hypercubic) lattice. 





Figure 1. Snapshot of the coarsening of ordered domains in the 2D Glauber-Ising model, after a quench 
to T = 1.5 < Tc from a totally disordered state and at times t — 25 (left) and t = 275 (right) after the 
quench. The black and white colours indicate the value of the Ising spins. 

A typical example of the microscopic evolution of such a system is shown in figure ^ Quite 
rapidly, ordered domains are formed which slowly move and grow. Empirically, it is found that 
the typical size of the domains scales with the time t as L(t) ~ t^^^, where z is the dynami- 
cal exponent. It is known that for dynamical rules chosen such that there is no macroscopic 
conservation law, z = 2 for quenches to T < Tc- The slow temporal evolution of macroscopic 
observables we are interested in results from the slow motion of the domain walls. It is common 
to consider the coarse-grained order parameter (e.g. mean magnetization for magnetic systems) 
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(f){t, r) and one tries to capture its time-evolution through a stochastic Langevin equation 

50 dn[^] 



dt 



+ r/ 



(1) 



where the gaussian white noise t] = r](t, r) is characterized by its first two moments {r]{t, r)) = 
and {i]{t, r)ri{t' , r')) = 2T5{t — t')5{r — r') and 7i is the Ginzburg-Landau functional. For systems 
undergoing a conventional second-order phase-transition, one expects that qualitatively 



n[4>\ = <p^<p + v[4>\ ; v[4>\ 



if T > Tc 
if T <Tr 



(2) 



where c^o = 4'o{T) are the two equilibrium values of (j) and A is the spatial Laplacian. Physically, 
systems with T > and T < are very different, since in the first case, there is a single 
ground-state (where Tilcj)] = min!) while there are two distinct ground states in the second case. 
Therefore, for T > the system will rapidly relax towards it single ground-state and no ageing 
occurs. On the other hand, if T < Tc it will depend on the microscopic environment of each 
spin variable to which of the two possible local ground-states the system will evolve locally. 
The competition between these distinct states then leads to ageing phenomena. 

The noisy Langevin equation can be turned into an equivalent field theory through the 
Martin-Siggia-Rose formalism, see e.g. [Sj. Schematically, it may be represented through the 
effective action 



s[<p,4>] 



dt dr 



9^ + 



5n\ 



where (p is the so-called response field conjugate to 
motion take the form 



(3) 



From this, the classical equations of 



(4) 



but initial conditions must still be specified. For interacting fields, fluctuation effects are not 
taken into account by the equations (0J), although they are present in the action S and the 
associated path integral. It turns out, see for a review, that ageing is more fully revealed 
through the study of two-time correlators C(t, s) and response functions R{t, s) defined by 



Cit,s) = m,r)<l){s,r)) ; 



5h{s, r) 



{<P{t,r)(f>is,r)) 



(5) 



h=0 



where h{s, r) is the local magnetic field at time s and the position r. The last equation comes 
from Martin-Siggia-Rose theory. Furthermore, = 0. 



Definition 1. A statistical system described by a Langevin equation or an effective action 
(jSJ is said to undergo ageing, if C = C{t, s) or R = R{t, s) depend on both the observation time 
t and the waiting time s and not merely on the difference t = t — s. 

If the times t, s and t — s become large simultaneously (as compared to some microscopic 
time scale tmicro)i one usually finds, for T < Tc, the following scaling behaviour, see e.g. 



C{t,s)^M^^fc{t/s) ; fc{x) 
R{t,s)c^s-^-^fR{t/s) ; fnix) 



~ X 



-Xc/z 



X 
X 



oo 
oo 



(6) 



where Meq is the equilibrium magnetization and Xc,R are the autocorrelation |^ and autore- 
sponse exponents, respectively. If long-range initial correlations of the form Cini{r) ~ 
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|r|~ are used, where d is the number of space dimensions and a a free parameter, one has for 
/erromagnets the rigorous bound Ac > (d + a)/2 (27j (it need not hold for disordered systems 
Furthermore, if a < 0, the relationship Xc = Xr + a has been conjectured [23]) while 
Ac = Xr for a fully disordered initial state is generally accepted. Finally, the value of the expo- 
nent a has recently been shown by scaling arguments to depend on the equilibrium spin-spin 
correlator Ceq as follows. Systems of class 5 have short-ranged correlators Ceq(r') ~ e~l^l/'» and 
systems of class L have long-ranged correlators Ceq{r) ~ Then |16j 

_ J 1/z ; class S , . 

"~ \ {d-2 + 7])/z ; class L ^ ^ 

This concludes our review of those main properties of ageing systems which we shall need below. 



2 Beyond scale invariance 

Our central question about ageing is the following: Is there a general, model-independent way 
to predict the form of the scaling functions fc{x) and fR{x) as defined in © ? 

The possibility of an affirmative answer might be suggested by the following known facts: (i) 
Ageing phenomena show a dynamical scaling, i.e. they are scale- invariant, (ii) In equilibrium 
critical phenomena, the invariance of the theory under dilatations r i— > 6r may be extended 
to local scale or conformal transformations r i— > b{r)r (such that angles are kept fixed). It is 
well-known that conformal invariance is a very powerful principle in two-dimensional equilibrium 
critical phenomena, see e.g. J3|- More precisely, we inquire whether the global dynamical scale 
transformations t ^ b^t and r ^ hr can be extended, analogously to conformal invariance where 
z = 1, to space-time dependent dilatation factors b = b{t,r) in a meaningful way. 

Example 1. Let z = 2 and consider d space dimensions. The Schrddinger group Sch{d) is 
defined by ^ 



at + (3 TZr + vt + a 

t ^ , r ^ — ; cx5-l3-i = l 8) 

7t + 7t + 

where TZ G SO{d), a,v M.'^ and a,/3, 7,5 S ffi. It is well-known that Sch{d) is the maximal 
kinematic group of the free Schrodinger equation Sip = with S = 2m\dt — I2(J| (that is, it 
maps any solution Sip = to another solution). Additional examples of Sch{d) as a kinematic 
group include certain non-linear Schrodinger equations systems of Schrodinger equations 
[S] and the Euler equations of fluid dynamics We denote the Lie algebra of Sch{d) by scl)^. 
Speciflcally, scf)^ = {X±i^q,Y^i/2, Mq} with the non-vanishing commutation relations 

[Xn, Xn'] = {n- n')Xn+n' , [Xn, Y^] = - m) Yn+m , [Yi/2, ¥.1/2] = Mq (9) 

where n,n' G {±1,0} and m € {±1/2}. 

Example 2. For a more general dynamical exponent z^2,we construct infinitesimal generators 
of local scale transformations from the following requirements (for simplicity, set d = 1): 
(a) Transformations in time are t i— > (of -|- (3)/{'^t + 5) with a5 — f3^ = 1. (b) The generator 
for dilatations is Xq = —tdt — z~^rdr — x/z, where x is the scaling dimension of the fields 
on which the generators act. (c) Space-translation invariance is required, with generator —dr- 
Starting from these conditions, we can show by explicit construction that there exist generators 
Xn-, n G {±1, 0} and Ym, m = — l/z, 1 — 1/z, . . . such that 

[Xn, Xn'] = (n - n)Xn+n' , [Xn, Ym] = (~ ~ ^n+m (10) 
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For generic values of z, it is sufficient to specify the 'special' generator 

X, = -t^dt - Ntrdr - Nxt - ar^d^-' - ^r^ dl^"" ' - T^^^^^-^^/^r^ (11) 

explicitly, where we wrote z = 2/N and 5,/?, 7 are free constants. Further non-generic solutions 
exist for = 1 and N = 2 In particular, the generator (|llj) reproduces for z = 2 those of 
Sch{l). The condition [Xi,y/v/2] = is only satisfied if either (I) /3 = 7 = which we call type 
/or else (II) 5 = which we call type II. 

Definition 2. If a statistical system is invariant under the infinitesimal generators of either 
type I or type II it is said to be locally scale-invariant of type I or type II, respectively. 

Only the generators of type II are suitable for applications to ageing phenomena. 

Theorem 1. The generators Xn,Ym of type II form a kinematic symmetry of the differential 
equation Sip = where 

S = -0 + ^)dt + ^d^, (12) 
z^ 

For z = 2, we recover the d = \ free Schrodinger equation and its maximal kinematic Lie algebra 
scf)i. See for the precise definition of the commuting fractional derivatives d^. 

In order to be able to apply this kinematic symmetry to ageing phenomena, we must consider 
the subset of {X^,!^} where time translations (generated by X_i = —dt) are left out. It can 
be checked that the initial line t = is kept invariant. Ageing systems are not time-translation 
invariant and local scale invariance for them is meant to exclude the generator 

Theorem 2. Consider an ageing statistical system which is locally scale-invariant of type II. 
Then the autoresponse function is 

Rit,s) = ro(-] {t-sy'-'' ; t>s (13) 

where tq is a normalization constant. Furthermore, consider the following scaling form of the 
spatio-temporal response 



6h{s,0) 

Then ^(u) is a solution of the equation 

du + z{P + j) udl-' + 2z{2 - z^dl" ] ^{u) = Q (15) 



In the special case z = 2, we have 

R{t,s;r) = Rit,s)e^pi- — ^] (16) 

where Ai = (3 + ^ is a constant. 

Proof. jl4| I15j Consider first the autoresponse R = R{t,s) = {(l)4>), where (p,(l) have the scaling 
dimensions x^f,, x^, respectively. Local scale invariance means that XnR = YmR = 0, with n > 0. 
Because of spatial translation invariance and the commutators (jlUj) . it is sufficient to check that 
XqR = XiR = 0. The explicit form <\11^ produces two linear differential equations for R{t,s) 
which are readily solved. Comparison with the scaling forms © then establishes p3p. Eq. H15|) 
is proven similarly, see jl5j . Eq. H16|) had been found earlier ^ 
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In these explicit forms of R{t, s; r) it is also assumed that the system is rotation-invariant. If 
that is not the case, A4 is no longer uniform, but becomes direction-dependent [17| . 

Eq. ((T^ has been tested and confirmed in several spin systems undergoing ageing, notably 
the kinetic Ising model with Glauber dynamics in d = 2, 3 through intensive simulations |14 | llfi j 
and the exactly solvable kinetic spherical model in ci > 2 dimensions [U lllll^ . for T <Tc, the 
random walk 7^ and finally (up to logarithmic correction factors) the 2D XY model with T < Tc 

and the 2D critical voter model j^. Furthermore, eq. ()16|) has been numerically confirmed 
in the Glauber- Ising model, again for d = 2,3 and T <Tc JZl- However, small corrections with 
respect to ()13|) were found at T = Tc in a two-loop e-expansion ^ and in a recent self-consistent 
study FT9^ of the time-dependent Ginzburg-Landau equation at T = and which improves on 
the approximate Ohta-Jasnow-Kawasaki theory. We refer to the literature for details. 



3 On the Schrodinger group 

Having reviewed the important phenomenological result of local scale invariance as given in The- 
orem |51 we now discuss in more detail how a dynamical symmetry such as local scale invariance 
might come about. We shall do this here for the case z = 2 and therefore discuss the Schrodinger 
group and Schrodinger- invariant systems For simplicity, we often set d = 1. 

Under the action of an element g S Scii(l) of the Schrodinger group, a solution (j){t, r) of 
the free equation {2Aidt — d^)cj) = where = im is fixed, transforms projectively, viz. 
(j){t,r) I— > {Tg(j)){t,r) = fg{g~^{t,r))cl){g^^{t,r)) with a known ,20] companion function fg. 

Following (Tn|, we treat M as an additional variable and ask about the maximal kinematic 
group in this case jJH]- First introduce a new coordinate C and a new wave function -0 through 

<A(t,r) = ^ / dCe-''^^V(C,t,r) (17) 
v2vr Jk 

We denote time t as the zeroth coordinate and C as coordinate number —1. When working out 
the action of the generators oi sci)i on the function ip, it is easily seen that the projective phase 
factors can be absorbed into certain translations of the variable C EHl- Furthermore, the free 
Schrodinger equation becomes 

In order to find the maximal kinematic symmetry of this equation, we recall that the three- 
dimensional Klein-Gordon equation d^d'^^{^) = has the 3D conformal algebra conf3 = 
so(4, 1) = B2 as maximal kinematic symmetry. By making the following change of variables 

c = ^(Co + ie-i) , t = ^(-^o + ie-i) , ^ = ^^1 (19) 

and setting Tp{(,t,r) = ^'(^), the 3D Klein-Gordon equation reduces to (fT8|) . Therefore jJH] 

Theorem 3. For variable masses Ai, the maximal kinematic symmetry algebra of the free 
Schrodinger equation in d dimensions is isomorphic to the conformal algebra conf^^2 ^'^'^ ^''^^ 
has the inclusion of the complexified Lie algebras {5ci)^)c C (conf(^^_2)c- 

For d = 1, the Cartan subalgebra is spanned by the generators Xq and N := —tdt + C^^C- 
With respect to the six generators of the Schrodinger algebra scf)^, there are the four additional 
ones N,V-^-,V-,W . They are identified in figure |2^ with the roots in the root space of i?2- 

Consider the non-isomorphic parabolic subalgebras of S2. In figure |21 these correspond to 
convex subsets of roots. There are two maximal parabolic subalgebras, namely (i) sclji := scf)]^© 
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Figure 2. (a) Roots of the complex Lie algebra B2 and the identification of the generators of the 
complexified conformal Lie algebra (conf3)c 3 (scf)]^)c- The double circle in the center denotes the 
Cartan subalgebra. The generators belonging to the three non-isomorphic parabolic subalgebras ^H] are 
indicated by the full points, namely (b) sciji, (c) oge^ and (d) alti. 



CN and (ii) a(ti := o[ii eCA^ where o[ti := {Xq, Xi,Y±i/2, Mq, V+} and V+ = -2Cr9^ - 2trdt - 
(r^ + 2i(t)dr — 2xr, see also figure|3)d. The minimal parabolic subalgebra is oge^ := oge^ © CA^, 
where agei := {Xq, Xi,Y±i/2, Mq}, see figure|2t. We see that both alti and agti do not contain 
the generator X_i of time translations and they may therefore be considered candidates for a 
dynamical symmetry of ageing systems. 

In writing down above the Klein-Gordon equation, we had used units such that the 'speed 
of light' c = 1. It had been claimed in the literature that in a non-relativistic limit c — > 00, 
there were a group contraction conf^^^ ^ ^'^^d provided that an ill-defined ". . . transfer of the 
transformation of mass to the coordinates ..." is carried out 2 . However, going over first to 
the variable (" before letting c — > 00, we do not find a group contraction but rather the map 
conf3 olti / scf)i. 

We now discuss some physical consequences of Theorem |31 Consider the effective action 
S = S[ip{C,t,r)]. For a free field one recovers from the action 



S= /dCdtdr 



^. dip dtp / dtp 
dC dt \ dr 



+ Sini (20) 



where Sini describes the initial conditions, the Schrodinger equation (|18|) as equation of motion. 
For an infinitesimal coordinate transformation parameterized by Eu, v = —1,0,1, a theory is 
said to be local if its action S transforms as (the second integral is restricted to the line t = 0) 



6S= /dCdtdrr/a'^e^+ / dCdrVe,, (21) 

J J{t=0) 
where T^'^ is the energy-momentum tensor . 

Theorem 4. // the action S of a local theory is invariant under translations in C and r, scale- 
invariant with z = 2 and Galilei-invariant, then 5xiS = and S in invariant under the action 
ofaQti- 



Proof. [iSI The line i = is invariant under the action of age^. Furthermore, locality (|2H) 
yields several Ward identities. First, translation-invariance in ^ and r implies U"^ = = 0. 
Dilatation invariance gives Tq^ + ^Ti^ = and Galilei-invariance implies Tq^ — iTi^^ = 0. Now, 
for an infinitesimal special Schrodinger transformation, 

6x,S = -e [dCdtdr [{2To° + Ti^) t + {Tq^ - iTr^) r] + '-^ [ dCdrr^U-^ = 
J 2 J(t=o) 



because of the Ward identities derived above. 



Dynamical space-time symmetry for ageing systems 



7 



For the free- field action (|20j) . an energy-momentum tensor which satisfies these Ward identities 
can be explicitly constructed ^Hj. We stress that Galilei- invariance must be included as a 
hypothesis and one cannot expect to be able to derive it from weaker assumptions. On the other 
hand, time-translation invariance is not required. 

Finally, we reconsider the derivation of the spatio-temporal response for z = 2 jl8j . A 
standard result of conformal invariance (see ^HI ) gives together with Theorem Ol 



1 r ^ 



(V'l(Cl,tl,ri)V2(C2,t2,r2)) = (^(^1)^(^2)) = ^0 5xi,.2 y^+2T) ^^^^ 

where xi^2 are the scaling dimensions of the fields '01,2) V'o is a normalization constant and 
C = Ci — C2 and so on. Using the convention > 0, it turns out that eqs. H17I22() together 

imply that {(j)i(f>2) = 0.^ However, if we define a conjugate field (ff by formal complex conjugation 
in ((T7I) with > 0, then for xi > we find 

(0i(ti,ri)0^(t2,r2)) = ^1 dCidC2e-'^^^^+*^2^MV'i(Ci,ti,ri)V'2(C2,t2,r2)) 

= (t>o5x,,xJ{Mi-M2)Q{t)t-^'ei,vi-^'-\ (23) 



2 t , 

and we stress that the 0-function expresses the causality condition ti > t2 which is required for 
an interpretation of {(j)i(f)2) as response function. If we recall from Martin-Siggia-Rose theory 
that R{t,s;r) = {cf>{t,r)(j){s,0)) = ((/>(t, r)0*(s, 0)), where effectively (j) has a 'mass' > and 
<j)* has a 'mass' 7W* < 0, we have the identification 

4>*{s,r) = ^{s,r) (24) 

This result has also been confirmed for the three-point response functions |18j . 

We now apply this to the parabolic subalgebras. For a system invariant under age^, we find 

/. \(x2-xi)/2 / ■ 2\-(a;i+2^2)/2 

{MCuh,n)MC2,t2,r2)) = 4^o [yj r(-i+-2)/2 [c + ^jj (25) 

If we would consider the extension agci — > sci)i, then invariance under time translations gen- 
erated by X-i = —dt fixes xi = X2 and we simply recover the result ()22() coming form the 
invariance under conf3. On the other hand, the extension oge^ — > alti requires that the condi- 
tion V+(V'iV'2) = holds. It is easy to see from ()25p that this is the case. Therefore, for both 
age^ and olti we find, provided xi + X2 > 



.21 = ct^oKMi - M2)m ( -i i-(-i+-2)/2 ) (26) 



{xi+X2)/2 

t2j ' ' V 2 t 



and we have indeed rederived (fT^ for ageing systems with z = 2 in a model-independent way, 
including^^the causality condition. There is not yet a criterion which would allow to decide if 
aQti or alti, if any, is the dynamic symmetry of ageing systems with z = 2. 

In conclusion, the explicit form of the spatio-temporal two-time response function R[t, s; r) as 
given in ()16() is a consequence of the assumed Galilei-invariance of the ageing system. The high- 
precision numerical confirmation of this form in the kinetic Ising model with Glauber dynamics 
quenched to T < Tc for d = 2 and d = 3 dimensions is a strong indication in favour of 
Galilei-invariance in that model. However, even for T = the Langevin equation or the 
system are for general 'H[(/>] not Galilei/Schrodinger-invariant, see [SlI^. In view of small 
corrections to (fT^ suggested by two- loop field-theory calculations 1^ and of Theorem |1] on 
the other hand, it remains to be understood in what precise sense ageing systems might be said 
to be Galilei- invariant. 



^This holds true for T = 0. For the case < T < Tc, see 
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